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ABSTRACT 

Various types of weather communications are required to alert 
industries and the general public about the impending occurrence of tor- 
nados, hurricanes, snowstorms, floods, etc. A natural disaster warn- 
ing satellite system has been proposed for meeting the communications 
requirements of the National Oceanic and Atmospheric Administration. 
Message traffic for a communications satellite was simulated with a 
digital computer in order to determine the number of communications 
channels to meet system requirements. IPoisson inputs are used for 
arrivals and an exponential distribution is used for service. 

INTRODUCTION 

The National Oceanic and Atmospheric Administration and the 
National Aeronautics and Space Administration have been jointly investi- 
gating various technologies in order to develop conceptual communica- 
tions systems which meet requirements for a natural disaster warning 
system. The function of such a system would be to: 

(1) Route disaster warnings to the general public. 

(2) Provide disaster communications among national, regional and 
local weather service offices and affected areas. 

(3) Provide environmental information to the general public. 

(4) Provide a system for collecting decision information for warning 
to the public. 
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The natural disasters which would be monitored by the disaster warn- 
ing system include tornados, severe thunderstorms, flash floods, tsunami, 
earthquakes, hurricanes, forest fires, winter storms, air pollution, etc. 

The National Weather Service is organized to monitor and predict the 
weather locally, regionally and nationally. There are also national centers 
for particular types of weather, for example, the National Hurricane 
Center in Miami, Florida. The total number of offices and centers around 
the country is approximately 300. 

The joint investigations by NOAA and NASA include terrestrial and 
satellite communication systems. This report is confined to a satellite 
system only. The problem is to determine the number of communications 
channels required for a satellite system. The information required for 
such a decision is difficult to generate since historical records show only 
the number and size of communications from various parts of the country. 
The exact time of transmissions cannot be determined and so it is impos- 
sible to determine instantaneous flows of- message traffic thus precluding 
a deterministic analysis of any network. Because of the local and regional 
nature of many communications, no individual has an intuitive understanding 
of the total problem. 

As will be demonstrated, the problem may be formulated as a multi- 
server queueing system. Simulation is frequently used toanalyze unique 
queueing-type problems which defy direct analytical solution. This tech- 
nique often provides more information than an analytical model because it 
is possible to formulate stochastic. simulation models which reveal the sys- 
tem states during occurrences of events with small probabilities of happen- 
ing, but which the system must be capable of handling. Such is the case of 
the natural disaster warning system. If messages are required to wait in 
a queue, a tornado may occur before the warning can be disseminated to the 
public. It is imperative that such a system would have minimum waiting 
times in a queue. 

The simulation model discussed in this report was formulated to handle 
the local, regional and national disaster warning communications of NOAA. 
If a Disaster Warning System were developed, it would be designed as an 
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interface with the many offices and centers throughout the country. The 
system would be used only to provide warnings to the public in the most 
expedient manner and to collect information from data collection platforms 
which would be located throughout the nation. The system would operate as 
an adjunct to the weather service rather than as a replacement for any 
present operation. 

The data collection platforms would be designed to monitor the environ- 
ment, for example, river and stream levels. This information would be 
relayed to a central area for data collection and then processed by the weather 
service. The channel allocation for such a system may be determined ana- 
lytically and so will not be treated here. Communication channels required 
for data collection platforms and teletypes may also be added to those 
determined necessary for voice communication messages. 

The classical queueing theory equations are discussed in this report in 
order to provide a framework for the development of a model; the equations 
are used to determine the expected values of certain parameters. 

CLASSICAL QUEUEING EQUATIONS (REF. 1) 

One of the most commonly encountered phenomena in the physical 
world is the waiting line process. The process occurs whenever a demand 
exceeds the capacity to provide service. In order to solve the waiting line 
problem, it is necessary to perform a trade-off between the "costs’ 9 of 
providing the service and the "costs" of not providing the service. Normally 
the goal is to achieve an economic balance between the two "costs" involved. 
Queueing theory and simulation models do not solve the problem directly, 
but the two approaches do provide the information required for decision 
making by predicting various characteristics of the queueing process. 

In the usual formulations of the process, units are generated over time 
by an "input source". These units enter the system and join a "queue". At 
certain points in time, a member of the queue is selected for service by some 
rule called a "service discipline". The required service is then performed 
for the unit by the "service mechanism", and then the unit leaves the queueing 
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system. The process is depicted in sketch (a). 

Queuei ng s y stem 

I 



(a) 


The size of the input source may be either finite or infinite. Since 
the calculations are easier for the infinite case, this assumption is often 
made even though the actual size is some relatively large finite number. 
The statistical pattern by which calling units are generated over time must 
also be specified. Usually it is assumed that this distribution is Poisson. 
An equivalent assumption is that the interarrival times form an exponential 
distribution since the cumulative distribution of the Poisson is of the expo- 
nential form 1-e 

The service discipline refers to the order in which members of the 
queue are selected for service. In this study it was assumed that the ser- 
vice discipline is first-come-first-served. 

The service mechanism consists of one or more facilities, each of 
which contains one or more parallel service channels or servers. The 
time elapsed from the beginning of service to completion is referred to as 
the service time or holding time. The probability distribution of service 
must also be specified for a queueing model. Special cases of the gamma 
distribution, the exponential distribution and constant service times are 
frequently selected for the service mechanism. 
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Although many types of waiting line situations have been studied, queue- 
ing theory has been primarily concerned with one particular situation, 
namely, a single waiting line with one or more servers as seen in sketch (b). 



The following is a listing of the standard notation and terminology 
used in queueing theory: 


Line Length 
Queue Length 


E 


n 


n 


number of calling units in the queueing system 

number of calling units waiting for service 

line length minus number of units being served 

state in which there are n calling units in the queueing 
system 

probability that exactly n calling units are in the queue- 
ing system 

number of servers or parallel service channels in the 
queueing system 





X n = mean arrival rate (expected number of arrivals per unit 

time) of new calling units when n units are in the sys- 
tem 

l± = mean service rate (expected number of units completing 

service per unit time) when n units are in the system 

L = expected line length 

L = expected queue length 

4 

W .-..=••• expected waiting time in the system (includes service time) 

W = expected waiting time in the queue (excludes service time) 

Si 

A negligible function of At or zero order effect will be denoted o(At). 

Since interest usually lies in a steady-state processes, rather than 
initial or startup conditions, queueing theory deals primarily with processes 
which are assumed to have reached a steady state. In this case, when X R 
is a constant, X, then 

L = XW 


and 


L, = XW n 

If the mean service time is assumed to be a constant, 1/j a then 

W = W„ + i 
q U 

The term ’’birth" refers to the arrival of a new calling unit into the 
queueing system and "death” refers to the departure of a served unit. Three 
postulates form the basis of the birth-death process. 

I. Birth Postulate: Given that the system is in state E n at time t, 
the probability that exactly one birth will occur in the interval from t to 
(t + At) is 

X fi At + o(At) 

where X^ is a positive constant. 
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II. Death Postulate: Given that the system is in state E n at time 
t, the probability that exactly one death will occur during the interval 
from t to (t + At) is 

M n At + °( At ) 


III. Multiple Jump Postulate: Given that the system is in state E n 
at time t, the probability that the number of births and deaths combined 
will exceed one during the interval from t to (t + At) is o(At). 

From the postulates it can be stated that one of four mutually exclu- 
sive events must occur during the interval from t to (t + At): 

1. Exactly one birth and no deaths. 

2. Exactly one death and no births. 

3. Number of births and deaths combined > one. 

4. No births or deaths. 

The sum of the four probabilities must equal one. The probability of 
event 4 equals 1- sum of probabilities for events 1 to 3, which during the 
interval from t to (t + At) is equal to 

1 - *- n At - M n At + o(At) 

since the sum or difference of o(At) terms can be written as o(At). The 
probabilities of being in state E n at time t + At are developed from the 
possible states at time t and the events required to go from that state to 
the state E as follows: 


State at t 


E 


J n-1 

n+1 

? 


E 


n 


Events from t to (t + At) 
one birth 
one death 
multiple events 


Probability of Occurrence 

p n -i (» n _i At + o(At)) 

P n+1 W + l At + °< At » 
o(At) 

P n (1 - A n At - (i n ULt + o(At) 


none 
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It is shown in reference 1 (p. 293) that 
dP 

— - = X 1 P -, + u P 1 - (X + a ) P for n > o 
^ n-1 n-1 ^n+1 n+1 v n ^n' n 

When n = o X 1 = o and u = o, so that 
n-1 o ’ 


^0 

dt 


= ^ 1 P 1 




This provides a set of differential equations which, if they could be 
solved, would provide the values for P . Unfortunately, a convenient 
general solution is not available and so the equations are used to obtain 
solutions for certain special cases. 


The Pure Birth Process 

Assume that X = X and u = o for all n>o. In this situation no 
n n 

deaths occur and the mean arrival rate is constant. The differential 
equations for this process are: 0 


^0 

dt 


- XP 0 


— - = XP , - XP for n = 1, 2, . . . 
dt n-1 n 


If the system is in state E Q at time t = o, then the solution for the n = o 
case is 
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The general solution is 


/>.,\n -Xt 

p _ (^i) ^ 

n " n! 

This is the Poisson distribution with parameter Xt. The mean and variance 

are both equal to Xt and the mean arrival rate is X. 

Although the pure birth processes not very interesting by itself, it does 

form one component of the queueing process used in many models. One of 

the results of this solution leads to a property referred to previously. 

P = e~ At implies that the probability that no births will occur during the 
O Xt 

time interval from o to t is e’ . Thus, the probability that the first 
birth will occur in this time interval is (1 - e - ^). If the random variable 
T is the time of the first birth then the cumulative distribution function of 
T is 

F(t) = p{t < t} = 1 - e _Xt , t > o 
Therefore, the probability density function of T is 

f(t) = Xe' Xt , t>o 

dt 

which is an exponential distribution. 

This result verifies that the expected time between arrivals is 



The Pure Death Process 

Assume that X r = o ^for all n > o and that M n = 4 for n > 1. Also 
assume that the system is in state E M at t = o. The first assumption 
implies that births never occur, and so this is a pure death process with 
a constant service rate until the process terminates at state E Q . The 



10 


results are similar to the pure birth process except that this process is 
the opposite. The differential equations reduce to 


dP 

— “ = ^n+l " 4P n for n = °* 2 > • • • M - 1 

dt 


dP 


M 


dt 


= -JLtP 


M 


M-n is the number of events that have occurred in this proccess. The 
probability that no events have occurred by time t is 



-Hit 


The probability that M-n events have occur ed 

P = ® for n = 1, 2, . . . , M 

n (M-n)! 

The remaining possibility is that M events have occured, so that 

M 

This is a truncated Poisson distribution with a parameter jut. The mean 
service rate is ju until the process terminates. The distribution of elapsed 
time between events is an exponential distribution. 


Steady State Solution 

The steady state solution for P may be obtained either by solving for 
P n in the transient case and letting t — oo or by setting dP n /dt = o in the 
differential equations and then solving for P . Since an elementary general 
transient solution is not available for the birth-death process, the second 
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approach will be used and an assumption made that a steady -state solution 
exists, i.e., 

lim P n (t) - P„ 


t~°° 


and 


lim 
t— °°| 


[dP n (t) 


dt 


= 0 


For the differential equations, 

0 = \-l P n-1 + ^n+l P n+1 ' (X n + V P n for n > 0 

o = /i* P, - X P for n = 0 
'l 1 o o 

The equation for n = 0 yields 


p i P ° 
4 1 


When n > o each equation yields 


X uP-X 1 P 1 

r> _ n t) , n n n-1 n-1 

n+1 " F n + 




n+1 


^n+1 


Considering the numerator of the second term when n > 1, 


/i - X 1 P 1 = u 

n n n-1 n-l n 


n-1 


4 


P i + 
n-1 


M n-1 P n-1 " X n-2 P n-2 


n 




n 


X n-1 P n-1 ^n-1 P n-1 ” X n-2 P n-2 



For successively smaller values of n this procedure must yield 


ll P - X , P 
'n n n-1 n- 




o o 


From the solution to the n = o equation 


^1 




P 

o 


so that 


ju P - X 1 P , = 0 
^n n n-1 n-1 


Then 


P -- n ”l p 
n *n-l 




n 


n-1 




n 


M 


n-2 

n-2 


n-2 


A. 4 X r> • * • A. 
n-1 n-2 o 




n ^n-l 


li. 


o 


or 


P 


n 



for n = 1, 2, 


To determine P , it is known that 
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so that 


P o = 


n-1 

~ 7T\ 

i + Z 

n=l n 

TT 

i=l 1 


For this information 


and 


L-In 


n=o 


n 


L q - | s <- s > p „ 


The summations do have analytic solutions for special cases, one of which 
is the multiple server model with Poisson input and exponential service. 

No other types of output have been solved for the case when S > 1. The 
state probabilities for the Poisson input -exponential service will be used 
to approximate the state probabilities for the simulation model. The model 
assumes that arrivals occur according to a Poisson input with parameter 
A and that the service time has an exponential distribution with mean (1/ju). 
The mean service rate for the system is dependent on the state of the system 
E n « The mean service rate per busy server is p. Therefore, the overall 
service rate must be n/i provided that n < S. If n > S, so that all servers 
are busy, p n = Sp. This is a special case of the birth death process with 
A n = A and 

P n = n M if o < n < S 

= S if n > S 
M 

If A < Sp, the mean arrival rate is less than the maximum mean service 
rate so that 
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Since — < 1, 

Sju 


the limit of the series 


so that 


and 




if n > S 
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Let 


X 



Then 



-S)P n 




Since 


p < l the limit of 




r X 
o\(J- 


it (i - pV 
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Wq=M 

X 

w = Wq + - 

L = Lq + - 

I 1 


Weather Service Message Traffic and Distributions 

The data for the message traffic was provided by the National Oceanic 
and Atmospheric Administration’s Environmental Research Laboratories in 
Boulder, Colorado. The data were divided into three types of inputs in order 
to develop distributions which could be utilized in the simulation model: 
hurricanes reaching the east coast of the U. S. ; weather warnings; and river 
for casts and warnings. 

The number of hurricanes reaching the east coast of the United States 
per year is a random variable having a Poisson distribution with X = 1. 9 
(ref. 2). This information was used to develop a hurricane simulation for 
100 years. A multiplicative congruential uniformly distributed random 
number generator was used to develop random numbers (ref. 3). These 
numbers were then mapped to a cumulative Poisson distribution in order to 
obtain the Poisson events. The hurricane simulation was used to develop a 
’’worst case” as an input for the communication satellite simulation model. 

In the 40th year, two hurricanes reached the eastern part of the U. S. on 
July 13. On July 14, another hurricane reached the east coast. Finally, 
on October 3 of the 40th year, one more hurricane reached the east coast. 

Using data from Hurricane Camille which occurred from August 12-14, 
1969, a Poisson distribution was predicated for hurricane message traffic 
(for satellite simulator) with an estimated parameter of X = 0. 019 per 
minute during hurricanes. This assumption, if incorrect, will not affect the 
model appreciably because the traffic for a hurricane is very small relative 
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to the other two types of message traffic. The effect is to increase the 
satellite channel requirements only during the periods mentioned above. 

The assumption also causes the results to be more conservative since the 
occurrence of three simultaneous hurricanes is a very remote possibility. 

The weather warning data were provided for the 72 months from 
January 1966 to December 1972. The data included the categories: tor- 
nadoes and severe storms; hurricanes; small craft and gales; forecasts 
for inland lakes; winter storm warnings; and other. 

A Poisson distribution was also predicated for the weather warnings. 

A Chi-square test was performed to determine the goodness of fit for a 
Poisson distribution with a = 0. 05. The 72 months of data yielded a 
parameter estimate of X = 0. 1454 per minute. In order to work with inte- 
gral data, the test was performed on the expected number of messages per 

hour which yielded an estimate of X = 8. 5. The experimental value for the 

2 

Chi-square statistic was 18.1. The value of Xq q 5 with 11 degrees of 
freedom was 19. 675 so that the hypothesis of a Poisson distribution for 
the weather warnings could not be rejected. 

A time series analysis was performed on the weather warning data in 
order to determine trends and seasonal variations. The data are shown in 
figure 1. The trend was recovered by using linear regression. If x is 
the number of years from 1965, and y is the average number of messages 
per month for the year x, then the expression 

y = 632 x + 4163 

may be used to estimate the value of the expected number of messages per 
month for a given year. * The correlation coefficient of the regression was 
r = 0. 96. Table I shows the seasonal variation in percentage of deviation 
from trend and figure 2 is a graph of the irregular variations in percentage 
of deviation from trend. 

The trend shows that the average number of monthly messages is in- 


1 


x = 0, 1, 2, . . . from base year 1966. 
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creasing at the rate of 632 per year. Therefore, the Poisson parameter 
should be increased in order to allow for a larger number of messages per 
month. It was not determined if there was actually more storms or whether 
there is a tendency to saturate the communications facilities, but the latter 
seems more likely. The Poisson parameter used for the simulation was 
based on the trend value for 1972 which yielded a value of A = 0. 1923 mes- 
sages per minute. 

The river forecast and warning data are treated in the same manner 
as the weather warning data. Data were obtained for the sixty months from 
January 1967 to December 1971. A Chi-square test was performed to 
determine the goodness of fit for a Poisson distribution with a = 0.05. The 
data yielded an estimate of X = 0. 5167 messages per minute. This was 
converted to 31 messages per hour. The hypothesis of a Poisson distribution 
could not be rejected at the oi = 0. 05 level. 

A time series analysis was performed on the river forecast data to 

determine the trend and seasonal variations in the same manner as was 

done for the warning data. Using the same notation as previously the aver- 

o 

age number of messages per month for the year x is given by 

y = 2667 x + 15665 

The correlation coefficient for this regression was r = 0.94. Table II 
shows the seasonal variation in percentage of deviation from trend and fig- 
ures 3 and 4 are graphs of the trend and irregular variations. 

The trend shows that the average number of messages per month is 
increasing at the rate of 2667 per year. The Poisson parameter was adjusted 
to allow for a larger number of messages per month based on the year 1972 
(X = 0. 7224 messages per minute). 


J x = 0, 1, 2, 


from base year 1967. 
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Message Processing Times 

A classification was made of 21 different types of weather service 
warnings and the average word length was provided by NOAA’s Environ- 
mental Research Lab in Boulder . The average length of all 21 types was 
136 words which also approximates the average speaking rate per minute. 

No data were given on the frequencies of the 21 message types but the 
average word length of each type was given. 

Since the parallel-channel queueing equations require exponential 
service, this distribution was selected arbitrarily. The average process- 
ing time equals approximately one minute assuming a speaking of 137 words 
per minute. It seems plausible that the majority of messages would require 
1 or 2 minutes to transmit, but that occasionally, messages would be on the 
order of 5 to 6 minutes. The exponential distribution allows for this possi- 
bility. If the parameter (i = 1 is used for the distribution, then the cumu- 
lative distribution of 1-e”^ where t is the processing time in minutes is 


Minutes 

Cumulative 

Delta 


probability 

probability 

1 

0. 632 

0. 632 

2 

.865 

.233 

3 

.950 

.085 

4 

.982 

.032 

5 

.993 

.011 

6 

.998 

.005 

7 

.999 

.001 

8 

1.000 

.001 


The delta probabilities may be interpreted to mean that 63. 2 percent of 
all messages will Mve a processing time of 1 minute; 23. 3 percent have 
times of 2 minutes; 8. 5 percent have times of 3 minutes, etc. Only inte- 
gral values were used for processing times to allow the computer program 
to perform most operations in integer arithmetic. 
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The Simulation Model and Computer Program 


As stated previously, the simulation model was developed to utilize 
Poisson input and an exponential distribution for service. The computer 
program utilized integer data when possible to minimize the CPU time. 

The queueing process input consisted of three types of message traffic: 
warning messages; river forecasts; disaster communications during and 
after hurricanes. The Poisson parameters used for these inputs were: 


Warning messages 
River forecasts 
Disaster communication 
for hurricanes 


X = 0. 1923 
X = 0.7224 

X = 0.057 from July 13-21 
= 0. 019 from Oct. 3-11 
= 0 Otherwise 


The exponential service parameter was the same for all three inputs 
(ju = 1.008). The program organization consisted of a main routine and 
12 subroutines. The source program names are: 

Main Routine NQAA - Serves as an executive routine and initializes 
some parameters. Prompts user for satellite channel capacity and a seed 
for the random number generator. Also contains a report generator. 

Subroutine MACHST - Sorting routine which determines the soonest 
available channel and then allocates that channel for use. 

Subroutine FILL - Routine which calls the message distribution and 
service routines and converts each non-zero event into a message queue 
for one week in increments of one minute. 

Subroutines NORDIS, RIVDIS, HURD IS - Routines which set Poisson 
parameters for each type of message. Each calls Poisson generator and 
then converts Poisson variable to an integral number of messages, (0, 1, 
2, etc.). These integral events are then returned to subroutine FILL. 

Subroutine MPROC - Routine which sets the parameter i± and calls 
the exponential distribution subroutine to obtain a service time. 

Subroutine GSERV - Routine which updates channel times and accumu- 
lates idle channel times and waiting times for messages in the queue. 
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Subroutine AVTIM - Routine which calculates average time and number 
of messages in the system. 

Subroutine AVUTIL - Routine which calculates average fractional channel 
utilization and the average time spent in the queue. 

Subroutine POISS - Routine which converts a uniformly distributed 
random number to a Poisson distributed random number. 

Subroutine EDIST - Routine which converts a uniformly distributed 
random number toan exponentially distributed random number. 

Subroutine Rand - Routine which generates uniformly distributed random 
numbers between zero and one using a multiplicative congruential technique. 

Using the convention that a given level may call only one subroutine at 
the next lower level and that control is always returned to the calling sub- 
routine, the flow of the program is depicted in sketch (c). 

A copy of the program appears in appendix A. Sample outputs are given 
in appendix B. 
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RESULTS OF SIMULATION AND CONCLUSIONS 


The simulation program was used to simulate one week for channel 
numbers ranging from 1 to 20. The results are shown in table III and the 
utilization factors are plotted in figure 5. 

Although the parameter used for the exponential distribution was 1.008, 
the average message processing time for all runs asymptotically approaches 
1. 6 because processing times less than 1 minute were not considered. The 
effect of this restriction was a reduction in ju. to 0. 625. Thus the data from 
the simulation runs is somewhat conservative. 


One of the essential requirements of the Natural Disaster Warning Sys- 
tem is that there be no delay in the transmission of warning messages. From 
the data in table III, this requirement means that the number of channels must 
be greater than eight if the average processing time is 1.6 minutes or more. 

The queueing equations were used to analyze the sensitivity of the model 
to changes in the parameter ju. The probability of being in state zero was 
calculated for channels numbering from 3 to 20. Using P Q , the probability 
of being in state (S + 1) was calculated for each number of channels from 


3 to 20. This probability Pg^ 


is the probability of a message transmission 


being delayed. Table IV shows the probabilities for X = 0.9717 per minute 


and \i = 1.008 per minute. The value X = 0. 9717 occurs only during the 


period of 3 simultaneous hurricanes. Table V shows the probabilities P Q 
and Pg + ^ for /i = 0. 625 per minute or a service time of approximately 
1.6 minutes. 


Although it is somewhat unrealistic to even consider such probabilities 
as 0. 0000001, the concept may be employed to mean an almost virtual 
certitude that the event will not occur in practice. To ensure that the sat- 
ellite system would never reach state (S + 1) the arbitrary criterion was 


established that Pg^ 


0.0000001 would determine the number of channels 


sufficient to meet the no-delay requirement. 


From tables IV and V it can be seen that S = 9 is sufficient for a service 


time which averages approximately one minute and S = 11 is sufficient for 
!± = 0. 625 or a service time which averages 1. 6 minutes. 
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The probabilities P Q and Pg + ^ were also calculated for average 
service times of 2 and 3 minutes. The resulting estimates for S were 
12 and 14, respectively. 

As a verification of the model, there was no statistically significant 
difference between the calculated Pg + ^ and the number of delays occurring 
for X = 0. 9717 and p. = 0. 625 (service time = 1.6 minutes) for S = 3 to 
S = 8 (table III). 

On the basis of the data used to establish the model a selection of 


S = 10 channels would offer a number sufficient to meet the requirements 
with a considerable safety margin. If such a choice were made table VI 
demonstrates the effects of power degradation on the accessibility of the 
satellite. 


The information in table VI may be used to conclude that if 10 channels 
were selected, the satellite could operate and be used effectively even with 
a 50 percent degradation in power or transmission capability since delays 
would be expected to occur at the average rate of 6 per 10 000 messages 
transmitted. Moreover the maximum delay would probably not exceed 
1 minute. 
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APPENDIX A 
COMPUTER PROGRAM 

The computer program was written in FORTRAN IV and executed on 
an IBM 360/67. The operating system TSS (Time Sharing System) allows 
terminal type interactive processing and so the program was written to 
be executed in a conversational mode. 
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0000100 C MAIN ROUTINE FOR THE COMMUNICATIONS SATELLITE 

0000200 C 

0000300 C SIMULATOR FOR THE DISASTER WARNING SYSTEM. 

0000400 C 
0000500 C 

0000600 DIMENSION I CHAN( 200 ) , I D LE( 200 ) , I TYPJ( 3 ) 

0000700 INTEGER HI 

0000800 INTEGERS IJNO(30000), I J I ND( 30000 ) , IMCHDX(30000), I WA I T( 30000), - 

0000900 1IPROC(30000), IQUE(llOOO) 

0001000 DATA I TYPJ/ ' HURR ' , ' WARN ' , 1 R I V . 1 / , I CHAN/ 2 00*1/ 

0091100 C 

0001200 WRITE(6,1032) 

0001300 1032 FORMAT ( ' ',T2,'IF WEEK = 1,HIT RETURN; OTHERW I SE TYPE 111') 

0001400 C 
0001500 C 

0001600 C INITIALIZE TIME PARAMETERS 

0001700 READ( 5 , 1001 ) I RND 

0001800 I F ( I RND . ME . 0 ) GO TO 50 

0001900 C 

0002000 K=0 

0002100 I M I N I T = 0 

0002200 IV/K“0 

0002300 C 

0002400 C 

0002500 C PROMPT FOR ENTRY OF NUMBER OF CHANNELS FOR SATELLITE 

0002600 C 
0002700 C 

0002800 WRITE (6,1000) 

0002900 READ (5,1001) NOCHAN 

0003000 C 
0003100 C 

0003200 C PROMPT FOR ENTRY OF A SEED TO START THE RANDOM 

0003300 C NUMBER GENERATOR. 

0003400 C 
0003500 C 

0003600 WRITE (6,1002) 

0003700 READ (5,1001) IGESS 

0003800 READ (5,1001) ISKIP 

0003900 C 

0004000 C EVENTS WILL BE GENERATED TO SIMULATE ARRIVALS FOR 

0004100 C 60 MINUTES PER HOUR, 24 HOURS PER DAY, FOR 7 DAYS. 

0004200 C THIS INFORMATION WILL THEN BE USED TO FORM A QUEUE 

0004300 C WHICH IS THEN PROCESSED. AFTER PROCESS I NG, SEVEN 

0004400 C MORE DAYS OF INFORMATION ARE GENERATED AND 

0004500 C PROCESSED. THIS PROCEDURE IS CONTINUED UNTIL 

0004600 C A YEAR HAS ELAPSED IN THE SIMULATION. 

0004700 C ,) 

0004800 C 

0004900 GO TO 100 

0005000 C 

0005100 50 READ(9, 1030) K, I M I N I T, I WK, NOCHAN, IGESS, I SK I P, Z, KKJ1, KKJ2 , KJ1, K J2 

0005200 READ (9,1031) ( I CHAN( I JK) , I JK = 1, 200 ) 

0005300 CALL RAND(Z, IGESS, KKJ1, KKJ2, KJ1, KJ2 > > 

0005400 C 
0005500 C 




0005600 C 
0005700 C 
0005800 C 
0005900 C 
0006000 C 
0006100 100 
0006200 
0006300 C 
0006400 C 
0006500 C 
0006600 C 
0006700 C 
0006800 
0006000 
0007000 
0007100 200 
0007200 C 
0007300 C 
0007400 C 
0007500 C 
0007600 C 
0007700 C 
0007800 
0007900 C 
0008000 C 
0008100 C 
0008200 
0003300 
0008400 
0008500 
0008600 
0008700 C 
0008800 600 
0008900 
0009000 
0009100 
0009200 
0009300 
0009400 
0009500 
0009600 
0009700 
0009800 
0009900 
0010000 650 
0010100 
0010200 
0010300 
0010400 
0010500 
0010600 
0010700 
0010800 C 
0010900 C 

0011000 c 


SUBROUTINE FILL GENERATES THE EVENTS AND TRANSFORMS 
NONZERO EVENTS INTO MESSAGES FOR A QUEUE. 


CALL FI LL(I JNO, I JIND, I PROC, K, IMIMIT, I JOB, I GESS, KK Jl, KK J2 , K Jl, K J2 ) 
K = 1 


SUBROUTINE GSERV PROCESSES THE QUEUE. 


DO 200 1=1, I JOB 

CALL GSERV( I , I JNO,MIN, I CHAN, NOCHAN, I MCHDX , I J IND, IDLE, IWAIT,- 
1 1 PROC, I QUE) 

CONTINUE 


SIMULATION FOR 1 WEEK COMPLETED. BEGIN PROCESSING 
OF QUEUE FOR PRINTING. 


PRINT HEADING ROUTINE 

WRITE (7,1004) 

WRITE (7,1005) 

WRITE (7,1006) 

WRITE (7,1007) 

WRITE (7,1008) NOCHAN 

I WK = I WK+ 1 
WRITE (7,1009) 

WRITE (7,1010) I WK 

MAXWT-0 

NOMAX =0 

NOARRV=0 

NOF I N=0 

DO 650 KK=1, I JOB 

IF (IJIND(KK).LE. 10080) NOARRV-NOARRV+ 1 
IFIN=IJIND(KK)+I WA IT(KK)+IPROC(KK) 

IF (IFIN.LT. 10080) NOF I N-NOF I N+l 

IF (( IWAIT(KK).GT.MAXWT).AND.(IJIND(KK).LE. 10080)) MAXWT = I WA I T( K K ) 

CONTINUE 

WRITE (7,1011) NOARRV 
WRITE (7,1012) NOF I II 
PDM I NS= 10080 
LO=l 

HI =10080 

CALL AVTI M( IJIND,LO,HI, NOARRV, ARRV, SUM, I JOB, IQUE,PDMI.NS,SUMQUE, IWAI T, I PROC) 
CALL AVUTI L( IDLE, SUM I DL, NOCHAN, SUMWT, I J IND, IWA IT, ARRV, H I , LO, I JOB, PDM I NS) 


AVTIM AND AVUTIL ARE USED TO CALCULATE THE AVERAGE 



0011100 

0011200 
0011300 
0011400 
0011500 
00 11C 00 
0011700 
0011800 
0011300 
0012000 
0012100 
0012200 
0012300 
0012400 
0012500 
0012GO0 
0012700 
0012800 
0012900 
0013000 
0013100 
0013200 
0013300 
0013400 
0013500 
0013600 
0013700 
0013800 
0013900 
0014000 
0014100 
0014200 
0014300 
0014400 
0014500 
0014600 
0014700 
0014800 
0014900 
0015000 
0015100 
0015200 
0015300 
0015400 
0015500 
0015600 
0015700 
0015800 
0015900 
0016000 
001G100 
0016200 
0016300 
0016400 
0016500 
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C TIME II! THE SYSTEM, AVERAGE HUMBER IN THE SYSTEM, 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


690 


C 

C 

c 

c 

c 

700 


FRACTIONAL CHANNEL UTILIZATION AND THE AVERAGE TIME 
IN THE QUEUE. USING THIS INFORMATION, THE AVERAGE 
PROCESSING TIME CAN BE DETERMINED BY SUBTRACTING THE 
AVERAGE TIME IN THE QUEUE FROM THE AVERAGE TIME IN 
THE SYSTEM. 


SUM = AVERAGE TIME IN THE SYSTEM 

SUMQUE= AVERAGE NUMBER IN THE SYSTEM 

SUM I DL= AVERAGE FRACTIONAL CHANNEL UTILIZATION 

SUMWT = AVERAGE TIME IN THE QUEUE 

AVPROC= SUM'SUMWT = AVERAGE PROCESSING TIME 


AVPROC=SUM-SUMWT 
WRITE (7,1024) SUM 
WRITE (7,1025) SUMQUE 
WRITE (7,1026) SUMIDL 
WRITE (7,1027) SUMWT 
WRITE ( 7, 1028) AVPROC 


DETERMINE MAXIMUM WAITING TIME AMD FREQUENCY 


DO 690 KK = 1, I JOB 

IF ( IWAIT(KK) .LT.MAXWT) GO TO 690 

NOMAX*NOMAX+l 

CONTINUE 

WRITE (7,1029) MAXWT, NOMAX 
IF (ISKIP.EQ.l) GO TO 810 
I PAGE=1 

WRITE (7,1014) IWK, I PAGE 
WRITE (7,1016) 

WRITE (7,1017) 

WRITE (7,1018) 

WRITE (7,1016) 

IF (ISKIP.EQ.l) GO TO 810 
DO 800 KK = 1, I JOB 

IF (IJIND(KK).GT. 10080) GO TO 800 
I PG=M0D( KK, 55 ) 

IF (IPG-.NE.O) GO TO 700 
I PAGE=I PAGE+1 
WRITE (7,1014) IWK, I PAGE 
WRITE (7,1016) . 

WRITE (7,1017) 

WRITE (7,1018) 

WRITE (7,1016) 


CONVERT ARRIVAL TIME TO DAY-HR-MIN FORMAT 


I ART = I J I ND( KK) 

I REM=MOD( I ART, 1440) 
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0016600 
0016700 
0016800 
0016000 
0017000 
0017100 
0017200 
0017300 C 
0017400 C 
0017500 C 
GC17GC0 C 
0017700 C 
0017800 
0017900 
0018000 
0018100 
C01S200 
0018300 
C018400 
0018500 
0018600 
0018700 
0018800 C 
0018000 C 

nniocoo c 

0019100 C 
0019 200 C 
0010300 
0019400 
0010500 800 
0010000 C 
0010700 C 
0019800 C 
0019900 C 
0020000 C 
0020100 C 
0020200 C 
0020300 C 
0020400 810 
0020500 
GC20COO 
QC2070Q 
0020800 
0020900 825 
0021000 
G0211C0 
0021200 
0021300 
0021400 
0021500 900 
0021000 C 
0021700 C 
0021800 
0021000 C 
0022000 C 


IF (1REM.EQ.0) 1 DAY1= l ART/1440 
IF (IREM.NE.O) IDAY1=I ART/1440 + 1 
I REM= I ART- ( I DAY1-1 ) *1440 
I R EM 1 =MOD ( I R EM , 6 0 ) 

IF (IREM1.EQ.0) I HR1 = I REM/ 6 0 
IF (IREM1.NE.0) I HR1 = I REM/60 + 1 

= I ART-( I DAY 1-1 1*1440 - (1HR1-1)*60 


COI1 VERT FINISH TIME TO SAME FORMAT 


I F I M = I JIND(KK) + IWAIT(KK) + IPROC(KK) 
I REM=MOD( I F I M, 1440) 

IF (IREM.EQ.O) I DA2 = 1 F I M/1440 
IF (IREM.NE.O) I DAY 2 = 1 FI fJ/1440 + 1 
1 REM= I F I N- ( I DAY2-1 ) * 1440 
I R EM 1 =MOD ( I !i EM ,60) 

IF (IREM1.EQ.0) I HR2 = I REM/CO 
IF (IREM1.NE.0) I HR2 = I REM/GO + 1 
HIM2=I FIN-( IDAY2-l)*1440-( IHR2-1)»60 
KJTYP=I JNO(KK) 


PRINT MESSAGE LOG 


WRITE (7.1023) KK, IDAY1, I HR1.MI Nl, I DAY2, I HR2 , M I N2 , I TYPJ ( K JTYP ) , 
llMCHDX(KK), IPROC(KK), IWAIT(KK) 

CONTINUE 


IF WEEK IS 52, THE PROGRAM 13 FINISHED; 
OTHERWISE ALL TABLES MUST BE CLEARED FOR 
THE NEXT WEEK. 


CONTINUE 

IF (IWK.EQ.52) GO TO 1500 
J = Q 

00 300 1=1, I JOB 

I JUO( 1 ) =0 
I JIMD( I ) =0 
I I1CHDX ( I ) =0 
I WA I T ( I ) =0 
I PROC( I ) =0 

CONTINUE 


CLEAR IQUE AND UPDATE I CHAN. 
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0022100 C 

C022200 DO 910 1=1,10080 

0022300 910 IO.UE( I ) =0 
0022400 DO 940 1=1,200 

0022500 I DLE( I ) =0 

0022600 940 I CHAN( I )=1 

0022700 CALL RAND(Z, I GESS, KKJ 1, KKJ2, KJ1, KJ2 ) 

0022800 935 CONTINUE 

0022900 WR I TEC 8, 1030) K, IMINIT, IWK, NOCHAN, IGESS, I SK I P, Z, KK Jl, KKJ 2, KJ1, K J2 

0023000 WR I TEC 8, 1031 ) C I CHANC I JK), I JK=1, 2 00) 

0023100 GO TO 1500 

0023200 C 

0023300 C FORMAT STATEMENTS 

0023400 C 
0023500 C 

0023600 1000 FORMAT (' ' , T2, ' ENTER NUMBER OF COMMUNICATION CHANNELS IN FORMAT 13') 
0023700 1001 FORMAT (13) 

0023800 1002 FORMAT (' 1 ,T2, 1 ENTER A RANDOM NUMBER BETWEEN 1 AND 999 IN FORMAT 13') 
0023900 1004 FORMAT ( ' 1 * , T2, * ') 

0024000 1005 FORMAT (' ' ,T56, 1 ********************' ) 

0024100 1006 FORMAT (' ',T56,'* COMM. SATELLITE *') 

0024200 1007 FORMAT (' ',TS6,'* SIMULATION *') 

0024300 1008 FORMAT ( ' 0 ' , T4 4, ' NO . OF CHANNELS 3 ' , T61, I 3 , T 7 5 , 'WEEKS SIMULATED 3 52') 
0024400 1009 FORMAT C-',T2,' ') 

0024500 1010 FORMAT (' ' , T25, ' WEEK ' , T30, I 2 ) 

0024600 1011 FORMAT (' \T52,'NO. OF ARRIVALS DURING PERIOD =',T84,I6> 

0024700 1012 FORMAT (' ',T52,'NO. OF MSGS COMPLETED THIS PD. =' ,T84, 1 6) 

0024800 1014 FORMAT ( ' 1 ' , T25, ' WEEK 3 ' , I 3, T100, ' PAGE 3 ' , I 4 ) 

0024900 1016 FORMAT (' \T25,' - - - 

0025000 1 — ') 


0025100 

1017 

FORMAT ( ' 

' , T 2 5 , ' | MSG 

1 ARRIVAL TIME 1 

FINISH TIME | MSG. 

0025200 


1 CHAN. 

1 PROCESS | 

WAIT | 

') 


0025300 

1018 

FORMAT (' 

' ,T25, ' | NO. 

| DAY 

HR. MIN. | 

DAY HR. MIN. | TYPE 

0025400 


1 ASSGN . 

| MINS. | 

MINS. | 

1 ') 


0025500 

1023 

FORMAT ( ' 

' ,T25, ' I ' , 15, 

T33, ' 1 ' , 1 

14, 14, 1 5,T50, 

' 1 ' , 1 4, 1 4, 1 5,T66, ' I ' , ' 


00 25 600 1T75, ' | ' , I 6,T8G, ' | ', I 6, T98, ' I ' , I 5, 3X, ' I ' ) 

0025700 1024 FORMAT (' ', T5 2, ' AVERAGE TIME IN SYSTEM 3 ',T84,F7.1) 

0025800 1025 FORMAT (' ' ,T52, 'AVERAGE NO. IN SYSTEM =',T84,F7.1) 

0025900 1026 FORMAT (' ' ,T52, 'AVERAGE FRACTIONAL CHAN. UT I L . 3 ' , T84, F8 . 2 ) 

0026000 1027 FORMAT (' ', T5 2, ' AVERAGE TIME IN QUEUE 3 ',T84,F7.1) 

0026100 1028 FORMAT (' ', T5 2, ' AVERAGE PROCESSING TIME =',T84,F7.1) 

0026200 1029 FORMAT ( ' 0 ' , T5 2, ' THE MAXIMUM DELAY OF', 15,' MINUTES OCCURRED ', I 5 , ' TIMES') 
0026300 1030 FORMATC 6 1 10, F15 . 12, 4 1 10 ) 

0026400 1031 FORMATC 200110) 

0026500 C 
0026600 C 
0026700 C 

0026800 1500 CONTINUE 
0026900 STOP 

0027000 END 
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0000100 

0C00200 

0000300 

0000400 

0000500 

0000600 

0000700 

COQOSOO 

0000900 


SUBROUTINE MACHS T (MIN, I MACH , MOMACH) 
DIMENSION I MACH( 200 ) 

Ml H=1 

DO 100 J=2, NOMACH 

IF ( I MACH (MIN) . LE . I MACH ( J ) ) GO TO 100 
Ml M=J 

100 CONTINUE 
RETURN 
END 
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0000100 

0000200 

0000300 

0000400 

0000500 

0000600 

0000700 

0000800 

0000900 

0001000 

0001100 

0001200 

0001300 

0001400 

0001500 

0001600 

0001700 

0001800 

9001300 

0002000 

0002100 

0002200 

0002300 

0002400 

0002500 

0002600 

0002700 

0002800 

0002900 

0003000 

0003100 

9003200 

0003300 

0003400 

0003500 

0003600 

0003700 

0003800 

0003900 

0004000 

0004100 

0004200 

0004300 

0004400 

0004500 

0004G00 

0004700 

0004800 

0004900 

0005000 

0005100 

0005200 

0005300 

0005400 

0005500 


SUBROUTINE FI LU I JNO, I J I NO, I PROC,K, 1M I N IT , I JOB, I GESS , KKJ1 , KK J2, K Jl, K J2 ) 

SUBROUTINE FILL GENERATES THE EVENTS AND TRANSFORMS . 

NONZERO EVENTS INTO MESSAGES FOR A QUEUE. 

INT£GER*2 I JNO, I J I ND, I PROC 

DIMENSION UNO (30000), IJIMD (30000), I PROC (30 000 ) 

IF(K.EQ.l) GO TO 200 
I JOB = 0 

DO 100 1=1,10080 
I M I N I T = I 

CALL N0RDIS(N0EVTS,IGESS,KKJ1,KKJ2,KJ1,KJ2) 

I F ( NOEVTS . EQ . 0) GO TO 25 
10 I JOB= I JOB+1 
I JNO( I JOB ) =2 
I J I ND( I JOB ) = I 

CALL MPROC(MSERV, I GESS , KK Jl, KKJ2, KJ1, KJ 2) 

I PROC ( I JOB ) =MSERV 

NOEVTS=NOEVTS-l 

I F ( NOEVTS . GT . 0 ) GO TO 10 

25 CONTINUE 

CALL Rl VD I S ( NOEVTS , I GESS, KK Jl, KK J2, KJ 1 , KJ 2 ) 

I F (NOEVTS . EQ . 0 ) GO TO 100 
30 I JOB= I JOB+1 

I JNO ( I JOB ) =3 

I J I ND( I JOB ) = I 

CALL MPROC(MSERV, I GESS , KKJ 1, KKJ2, KJ 1, KJ 2 ) 

I PROC( I JOB) =MSERV 
NOEVTS=NOEVTS-l 
I F ( NOEVTS .GT . 0) GO TO 30 
100 CONTINUE 
110 GO TO 400 
200 DO 300 1=1,10080 
IMI NIT-IMI NIT+1 

I F(( IMINIT.LT. 279360). OR. (IMINIT.GT. 293760)) GO TO 210 
CALL HURD I S ( NOEVTS , 1 M I N I T , I GESS, KKJ 1, KK J2 , K Jl, KJ2 ) 

I F ( NOEVTS . EQ. 0) GO TO 210 
205 I JOB= I JOB+ 1 
I JNO ( I JOB ) =1 
I J 1 NO ( I JOB) = I 

CALL MPROC(MSERV, I GESS, KKJ 1, KKJ2 , KJ 1, KJ 2 ) 

I PROC( I JOB) =MSERV 

NOEVTS=NOEVTS-l 

I F( NOEVTS . GT . 0) GO TO 205 

210 IFUIMINIT.LT. 397440). OR. (IMINIT.GT. 411840)) GO TO 220 
CALL HURD I S ( NOEVTS, IMINIT, I GESS, KK J 1, KK J2 , K Jl, KJ 2 ) 

• IF(NOEVTS.EQ.O) GO TO 220 
215 I JOB= I JOB+ 1 
I JNO( I JOB ) =1 
I J I ND ( I JOB ) = I 
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0005600 

0005700 

0005800 

0005900 

0006000 

0006100 

0006200 

0006300 

0006400 

0006500 

0006600 

0006700 

0006300 

0C06900 

0007000 

0007100 

0007200 

0007300 

0007400 

0007500 

0007600 

0007700 

0007800 

0007900 

0008000 

0008100 

0008200 

000830C 

0008400 

0008500 


CALL MPROCCMSERV, I GESS, KK J 1, KKJ 2 , KJ 1, KJ 2 ) 
IPROCC I JOB)=MSERV 
NO E VTS =NO E VTS - 1 
IF(NOEVTS.GT.O) GO TO 215 
220 CALL NORDIS(MOEVTS, I GESS, KKJ1, KKJ2, KJ1, KJ2 ) 
IF(NOEVTS.EQ.O) GO TO 230 
225 I JOB = I JOB + 1 
I JNOC I JOB) =2 
I J I ND( I JOB) = 1 

CALL MPROC(MSERV, l GESS, KK J 1, KKJ2, K.J1, KJ2 ) 

I PROC( I JOB ) =MS ERV 

NOEVTS=NOEVTS-l 

I F ( NOEVTS . GT . 0 ) GO TO 225 

230 CONTINUE 

CALL R I VO I S ( MO E VTS , I GESS, KKJ1, KKJ2, KJ1, KJ2 ) 
IF(MOEVTS.EQ.O) GO TO 250 
235 I JOB= I JOB+ 1 
I JNO ( I JOB ) =3 
I J I MD( I JOB ) = I 

CALL MPROC (MS ERV, I GESS, KKJ1, KKJ 2, KJ1, KJ2 ) 

I PROC ( I JOB ) =MS ERV 
NO E VTS = NOE VTS -1 
I F ( NO EVTS . GT . 0 ) GO TO 235 
250 CONTINUE 
300 CONTINUE . 

400 CONTINUE 
RETURN 
END 
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0000100 
0000200 C 
0000300 C 
0000400 C 
0000500 
0000600 C 
0000700 C 
0000800 C 

0000000 c 
0001000 
0001100 
0001200 
0001300 
0001400 
0001500 


SUBROUTINE R I VO I S ( NOEVTS, IGESS / KKJ1,KKJ2,KJ1,KJ2) 


THIS SUBROUTINE CALLS A POISSON DISTRIBUTED GENERATOR 
WHICH DEVELOPS THE NUMBER OF SIMULTANEOUS RIVER 
WARNING MESSAGES 

REAL LAMDA 
LAMDA=0. 7224 

CALL PO I SS ( LAMDA, NOEVTS , I GESS, KK Jl, KKJ 2, KJ1, KJ2 ) 

CONTINUE 

RETURN 

END 
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0000100 
0000200 C 
0000300 C 
0000400 C 
0000500 C 
0000600 C 
0000700 C 
0000800 C 
0000900 C 
0001000 
0001100 
0001200 
0001300 
0001400 
0001500 
0001600 
0001700 
0001800 


SUBROUTINE HURD I S ( NOEVTS , I M I M I T, I GESS, KK J 1, KKJ2, KJ1, KJ2 ) 


THIS SUBROUTINE CALLS A POISSON DISTRIBUTED GENERATOR WHICH 
DEVELOPS THE SIMULTANEOUS NUMBER OF HURRICANE WARNING 
MESSAGES 


REAL LAM DA 

I F C IMINIT.GT. 293760) GO TO 1 
LAMDA= 0.057 
GO TO 2 

1 LAMDA-0.019 

2 CALL POI SS( L AMD A, NOEVTS , 1 GESS, KK Jl, KKJ2 , KJ1, KJ2) 
CONTINUE 

RETURN 

END 
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0000100 
0000200 C 
0000300 C 
0000400 C 
0000500 C 
0000G00 C 
0000700 C 
0000800 C 
0000900 C 

0001000 c 

ononoo 

0001200 

0001300 

0 0 C 1 4 9 0 

0001500 

0001600 


SUBROUTINE NORD I S ( NOEVTS , I GESS, KK J 1, KK J2 , KJ 1, K J 2 ) 


THIS SUBROUTINE CALLS A POISSON DISTRIBUTED GENERATOR 
WHICH DEVELOPS THE NUMBER OF SIMULTANEOUS MESSAGES OR 
EVENTS FOR A PROCESS WHICH HAS A PROBABILITY DENSITY 
FUNCTION WHICH -IS DISTRIBUTED AS A POISSON DENSITY 
FUNCTION. 

REAL LAMDA 
LAMDA=0. 1923 

CALL PO I SS( LAMDA, NOEVTS, I GESS , KKJ 1, KKJ 2, KJ 1, KJ 2 ) 

CONTINUE 

RETURN 

END 
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0000100 
0000200 C 
0000300 C 

0 0 0 0 1( 0 0 c 

0000500 C 
0000600 C 
0000700 C 
0000800 C 
0000900 C 
0001000 
0001100 
0001200 
0001300 
00011(00 
0001500 


SUBROUTINE MPROCCMSERV, I GESS, KKJ 1, KKJ2 , KJ 1, KJ 2 ) 


THIS SUBROUTINE IS USED TO DEVELOP A PROCESSING TIME FOR 
MESSAGES. THE DISTRIBUTION TIME IS EXPONENT I AL AND 
BASED ON HISTORICAL VALUES FOR MESSAGE PROCESSING TIMES. 


REAL MU 
MU-1.008 

CALL EOIST (MU.MSERV, I GESS, KKJ1, KILJ2, KO 1, KJ2 ) 

CONTI NUE 

RETURN 

END 
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0000100 
0000200 C 
0000300 C 
0000400 C 
0000500 C 
0000600 
00CC700 
0000800 C 
0000900 
0001000 
0001100 
0001200 
0001300 
0001400 
0001500 
0001600 
0001700 
0001800 
0001900 
0002000 
0002100 
0002200 
0002300 1105 
0002400 1110 

0002500 
0002600 
0002700 
0002800 1115 
0002900 
0003000 
0003100 
0003200 
0003300 
0003400 
0003500 1120 
0003600 
0003700 


SUBROUT I HE GSERV( I , UNO, MIN, I CHAM, NOCHAN, IMCHDX, I JIMD, IDLE, IWAIT, IPROC, IQUE) 
SUBROUTINE GSERV PROCESSES THE QUEUE. 


DIMENSION I CHAN{ 200 ) , IDLE(200) 

I NTEGER*2 UNO, IMCHDX, I JIMD, (WAIT, IPROC, IQUE 

DIMENSION IJ NO (30000), I MCHDX( 30000 ) , IWA I T( 30000 ) 
DIMENSION I PROC( 30000 ) , I QUE( 11000 ) , I J I ND( 30000) 

CALL MACHSTWH, ICHAH, NOCHAN) 

I MCHDX ( I ) =M I N 

IF( ICHAN(MIN) .GT. IJI !.'D( I ) ) GO TO 1105 
IWA I T ( I ) =0 


IDLE(MIN)=IDLE(MIN)+IJIND( I ) - I CHAN(MIN) 
I CHAN (M I N ) = I J I NO ( I ) 

GO TO 1110 

I WA I T ( I ) = I CHAN (M I N ) - I J I NO ( I ) 

I F I N= I J I ND( I ) + 1 WAI T( I )♦ I PROC( I ) 

I F I N 1 ■ I F I N 
IFIN2-I JIHD( I ) 

DO 1115 IT I ND = I F I N2, IFIM1 
IQUE( I T 1 ND ) = I QUEC I T I ND ) + 1 
I CHAN (M I N) = ICHAN(MIN) + I PROC ( I ) 


CONTINUE 

RETURN 

END 



GO 0 0 1 0 C 
0000200 C 
0000300 C 
0000400 C 
0000500 C 
0000600 
0000700 
0000800 
0000900 
0001000 
9001100 
0001200 
0001300 
0001400 
0001500 100 
0001600 
0001700 
0001800 
C001900 200 
0002000 
0002100 
0002200 


SUBROUTINE AVT I M( I J I UD , LO , H I , NOARRV, ARR V, SUM, I JOB, I QUE, PDM I NS , SUMQUE, I WAIT, I PROC) 

THIS SUBROUTINE CALCULATES THE AVERAGE TIME IN THE SYSTEM 
AND THE AVERAGE NUMBER OF MESSAGES IM THE SYSTEM. 

DIMENSION I J t ND (30000), IQUE(llOOO), I WA I T( 30000 ) , I PROC (30 000) 

I NTEGER*2 I JIND, IQUE, IWAIT, IPROC 

INTEGER HI 

SUMQU£=0. 

SUM=0. 

DO 100 KKK=1, I JOB 
IF (I JIND(KKK).LT.LO) GO TO 100 
IF (IJIND(KKK).GT.HI) GO TO 100 
SUM=SUM+ IWAI T (KKK)+ I PROC(KKK) 

CONTINUE 
ARRV=NOARRV 
SUM=SUM/ARRV 
DO 200 I =LO, H I 
SUMQUE=SUMQUE+IQUE( I ) 

SUMQUE =SUMQUE/ PDM I NS 

RETURN 

END 
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0000100 

0000200 

0000300 

0000400 

0000500 

0000600 

0000700 

0000300 

0000200 

0001000 

0001100 

0001200 

0001300 

0001400 

0001500 

0001600 

0001700 

0001800 

0001900 

0002000 

0002100 

0002200 


SUBROUT I ME AVUTI L( I OLE, SUM I D L, NOMACH, SUMWT, I Jl fiD, I WA I T, ARRV, H I , LO, I JOB, PDHRS) 
C THIS SUBROUTINE IS USED TO CALCULATE THE FRACTION OF 

C TIME THE COMMUNICATION CHANNELS ARE USED AND THE AVERAGE . 

C WAITING TIME IN THE QUEUE. 

C 

I NTEGER* 2 I J I ND, 1 VI A I T 

DIMENSION I DLE ( 200 ) , I J I ND{ 3 0000 > , IWAI T ( 30000 ) 

INTEGER HI 
SUM I DL = 0 . 

SUMWT=0 . 

DO 100 KKK=1, NOMACH 
100 SUM I 0 L=SUM I DL+ I D LE ( KKK ) 

DMACH=NOMACH 

SUM I D L = ( PD HRS -SUM I DL/DMACH ) / PDHRS 
DO 200 I =1, 1 JOB 
IF (IJINDCI ) . LT. LO) GO TO 200 
IF (IJINDCI ) . GT.H I ) GO TO 200 
SUMWT =SUMWT + I WA I T ( I ) 

200 CONTINUE 

SUMWT =SUMWT /ARRV 

RETURN 

END 
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0000100 
0000200 C 
0000300 c 

ooooaoo c 

0000500 C 

0000600 c 

0000700 
0000800 
0000900 
0001000 
0001100 
0001200 
0001300 
0001400 100 
0001500 
0001600 
0001700 
0001800 
0001000 
0002000 
0002100 
0002200 
0002300 
0002400 200 
0002500 300 
0002600 
0002700 


SUBROUTINE PO I SS( LAMDA, NOEVTS, I GES S , KKJ1, KK J 2 , KJ 1 , K J 2 ) 

THIS SUBROUTINE MAPS A UNIFORMLY DISTRIBUTED RANDOM NUMBER ONTO 
A CUMULATIVE POISSON DISTRIBUTION IN ORDER TO OBTAIN A POISSON 
DISTRIBUTED RANDOM NUMBER. 

DIMENSION PROBC10) 

REAL NFACT, LAMDA 

NFACT-1.0 

PZERO = EXP ( -LAMDA) 

DO 100 11=1,10 
NFACT=NFACT*N 

PROB(N)=(LAMDA»»N)*£XP(- LAMDA) /NFACT 
CONTINUE 

CALL RAND( Z, I GESS, KK J 1, KK J 2, KJ1, KJ 2 ) 

HOEVTS=0 
Z=Z- PZERO 

IF (Z.LT.0.0) GO TO 300 
NOEVTS=NOEVTS+l 
DO 200 N = l, 10 
Z=Z-PROB(N) 

IF (Z.LT.0.0) GO TO 300 

NOEVTS=NOEVTS+1 

CONTINUE 

CONTINUE 

RETURN 

END 



42 


0000100 

0000200 C 

0000500 C 

0000400 C 

0000500 C 

0000600 C 

0000700 

0000800 

0000900 

0001000 

9001100 

0001200 

0001300 100 

0001400 200 

0001500 

0001600 

0001700 

0001800 

0001900 

0002000 300 

0002100 

0002200 


SUBROUTINE ED I ST ( MU,MSERV , I GESS , KK J 1 , KK J 2, K J1 , K J2 ) 

THIS SUBROUTINE MAPS A UNIFORMLY DISTRIBUTED RANDOM NUMBER ONTO 
A CUMULATIVE EXPONENTIAL DISTRIBUTION IN ORDER TO 'OBTA I N AN 
EXPONENTIALLY DISTRIBUTED RANDOM NUMBER. 

DIMENSION PROB( 150 ) 

REAL MU 
DATA 1/1/ 

IF (I . EQ. 0 ) GO TO 200 
I =0 

DO 100 N=l,50 
PROB(N)=1.0-EXP(-MU*N) 

CONTINUE 

MSERV=1 

CALL RAND( Z , I GESS, KKJ 1, KK J 2, KJ1, KJ2 ) 

DO 300 N=l,50 

IF (Z.LT.PROB(N)) GO TO 300 

MSERV-MSERV+1 

CONTINUE 

RETURN 

END 
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0000100 

0000200 

0000300 

0000400 

0000500 

0000600 

0000700 

0000800 

0000900 

0001000 

0001100 

0001200 

0001300 

0001400 

0001500 

0001600 

0001700 

0001800 

0001900 

0002000 

0002100 

0002200 

0002300 

0002400 


SUBROUTINE RAND( Z , I GESS, A, X, I , I SW) 

C 

C SUBROUTINE RAND GENERATES UNIFORMLY DISTRIBUTED RANDOM NUMBERS 

C 

INTEGER A, X 

M=2**20 

FM=M 

IF (I .EQ.l) GO TO 100 
1=1 

X=566387 

A=2**10+3 

100 X=MQO(A»X,M) 

FX = X 
Z=FX/FM 

IF (ISW.EQ.l) GO TO 300 
DO 200 K=l, I GESS 
X=MOD( A*X,M) 

FX=X 

Z=FX/FM 

200 CONTINUE 
ISW = 1 

300 CONTINUE 
RETURN 
END 
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APPENDIX B 

SAMPLE COMPUTER OUTPUTS 

The output from the computer program consists of statistics and a 
message log. The first example consists of statistics and one page of 
the message log for week one of a simulation of 4 communication channels 
The second example consists of only the statistics for week 43 of a Simula 
tion of 10 communication channels. The message log may be printed or 
suppressed at the users option. 




* COMM. SATELLITE * i 

* SIMULATION * : 

NO. OF CHANNELS 0 4 WEEKS^ SIMULATED 15 52 


WEEK 1 • • 

NO. OF ARRIVALS DURING PERIOD - 9345 

NO. OF MSGS COMPLETED THIS PD.= 9343 
AVERAGE TIME IN SYSTEM = 1.6 

AVERAGE NO; IN. SYSTEM *> 2.4 

AVERAGE FRACTIONAL CHAN. UTIL.= 0.37 

AVERAGE TIME IN QUEUE * 0.0 

AVERAGE PROCESSING TIME = 1.6 

THE MAXIMUM DELAY OF 4 MINUTES OCCURRED 




•;u » * 



1 TIMES 
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VJEEK= 1 PAGE= 1 


1 

1 

MSG 1 

NO. I 

ARRIVAL 
DAY HR. 

TIME 

M 1 N .. 

F I M 1 SH 
DAY HR. 

TIME 

MIN. 

MSG. 

TYPE 

1 

1 

CHAN. I 
ASSGN . I 

PROCESS 

MINS. 

1 WAIT 
1 MINS. 

1 

1 

1 | 

1 

1 

1 * 

1 

1 

2 

WARN 

1 

- 1 1 

1 

1 0 

1 

1 

2 1 

1 

1 

1 ’• 

1 

1 

2 

R 1 V . 


; 2 I 

1 

1 0 

1 


3 | 

1 

1 

1 

1 

1 

2 . 

. RIV. 

1 

, 3 1 

1 

1 0 

1 

1 

4 | 

1 

1 

2 

1 

1 

3 

RIV. 

1 

4 | 

1 

1 0 

1 

1 

5 I 

1 

1 

4 

1 

1 

7 

WARN 

1 

1 I 

3 

1 0 

t 

1 

0 1 

1 

1 

4 

1 

1 

" 6 

' RIV. 

I 

2 1 

2 

1 0 

1 

1 

7 | 

1 

1 

5 

1 

1 

6 

RIV. 

1 

3 1 

1 

1 0 

1 

1 

8 1 

1 

1 

8 

1 

1 

9 

RIV. 


4 1 

1 

1 0 

1 


9 I 

1 

1 

9 

1 

1 

11 

RIV. 


2 1 

2 

1 0 

1 

1 

10 1 

1 

1 

12 

1 

1 

13 

RIV. 

1 

3 i 

1 

1 0 

1 

1 

11 1 

1 

1 

13 

1 

1 

14 

WARN 

1 

1 I 

1 

1 0 

1 

1 

12 I 

1 

1 

13 

1 

1 

14 

RIV. 

1 

4 1 

1 

1 0 

1 

1 

13 ! 

1 

1 

14 

1 

1 

IS 

RIV. 

1 

2 1 

1 

1 0 

1 

1 

14 I 

1 

1 

15 

1 

1 

16 

WARN 

1 

3 I 

1 

1 0 

1 

1 

15 I 

1 

1 

15 

1 

1 

16 

RIV. 


1 I 

1 

1 0 

1 


16 | 

1 

1 

15 

1 

1 

" 16 

RIV. 


4 1 

1 

1 0 

1 

1 

17 I 

1 

1 

16 

1 

1 

18 

WARN 

1 

2 1 

2 

1 0 

1 

1 

18 I 

1 

1 

18 

1 

1 

19 

WARN 

1 

1 1 

1 

1 0 

1 

1 

19 I 

1 

1 

18 

1 

1 

19 

RIV. 

1 

3 I 

1 

1 0 

1 

1 

20 | 

1 

1 

19 

1 

1 

20 

RIV. 

' | 

4 1 

1 

1 0 

1 

1 

21 1 

1 

1 

20 

1 

1 

21 

RIV. 

1 

2 1 

1 

1 0 


1 

22 | 

1 

1 

22 

1 

1 

23 

RIV. 

1 

1 1 

1 

1 0 

1 


23 | 

1 

1 

24 

1 

1 

25 

RIV. 

1 

3 1 

1 

1 0 

1 

1 

24 | 

1 

1 

28 

1 

1 

29 

RIV. 

1 

4 I 

1 

1 0 

1 

1 

25 | 

1 

1 

31 

1 

1 

32 

RIV. 

1 

2 | 

1 

1 0 

1 

1 

26 | 

1 

1 

32 

1 

1 

34 

RIV. 

1 

1 1 

2 

1 0 

1 

1 

27 | 

1 

1 

35 

1 

1 

36 

RIV. 

1 

3 I 

1 

1 o 

1 

1 

28 I 

1 

1 

36 

1 

1 

39 

WARN 

1 

4 1 

3 

1 0 


1 

29 I 

1 

1 

37 

1 

1 

40 

RIV. 

1 

2 I 

3 

1 o 

1 

1 

30 1 

1 

1 

38 

' 1 

1 

39 

RIV. 


1 1 

1 

1 o 

! 

1 

31 1 

1 

1 

39 

1 

1 

42 

WARN 

1 

3 1 

3 

1 o 

1 

1 

32 | 

1 

1 

44 

1 

1 

48 

RIV. 

1 

1 1 

4 

1 o 

1 

I 

33 | 

1 

1 

46 

1 

1 

49 

RIV. 

1 

4 1 

3 

1 o 

1 

1 

34 I 

1 

1 

46 

1 

1 

47 

RIV. 

X 

2 1 

1 

1 o 


1 

35 | 

1 

1 

46 

1 

1 

4 7 

R 1 V . 


3 1 

1 

1 o 


1 

36 I 

1 

1 

47 

1 

1 

48 

RIV. 


2 1 

1 

1 0 


1 

37 1 

1 

1 

51 

1 

1 

54 

WARN 

1 

3 1 

3 

1 o 

1 

1 

38 | 

1 

1 

51 

1 

1 

54 

RIV. 

1 

1 I 

3 

1 0 

1 

1 

39 | 

1 

1 

52 

1 

1 

53 

RIV. 

1 

2 1 

1 

1 0 

1 

1 

40 | 

1 

1 

52 

1 

1 

53 

RIV. 


4 1 

1 

1 0 

1 

1 

41 1 

1 

1 

53 

1 

1 

54 

WARN 

1 

2 1 

1 

1 0 

1 

1 

42 | 

1 

1 

54 

1 

1 

57 

WARN 

1 

4 1 

3 

1 0 

1 

1 

43 | 

1 

1 

56 

1 

1 

57 

RIV. 

1 

1 1 

1 

1 0 

1 

1 

44 | 

1 

1 

57 

1 

1 

58 

RIV. 


2 1 

1 

1 ' 0 


1 

45 | 

1 

1 

58 

1 

1 

59 

WARN 

1 

3 1 

1 

1 o 


1 

46 | 

1 

1 

58 

1 

1 

59 

WARN 

1 

1 1 

1 

1 0 

1 

1 

47 | 

1 

1 

58 

1 

1 

59 

RIV. 

1 

4 1 

1 

1 o 

1 

1 

48 | 

1 

1 

59 

1 

1 

60 

RIV. 

1 

2 1 

r 

1 o 

1 

1 

49 | 

1 

1 

60 

1 

2 

4 

RIV. 

1 

1 I 

4 

1 o 

1 

1 

50 | 

1 

2 

4 

1 

2 

5 

RIV. 

1 

3 I 

1 

1 o 

1 

1 

51 1 

1 

2 

5 

1 

2 

9 

WARN 

1 

4 1 

4 

1 o 

1 

1 

52 1 

1 

2 

5 

1 

2 

7 

RIV. 

1 

2 ! 

2 

1 o 

1 

1 

53 1 

1 

2 

5 ' 

1 

2 

6 

RIV. 

1 

1 1 

1 

1 o 

1 

1 

54 | 

1 

2 

6 

1 

2 

7 

WARN 

1 

3 1 

1 

1 o 

1 



I 


******************** 

* COMM. SATELLITE ':. 

* , SIMULATION * 

NO. OF CHANNELS® 10 WEEKS SIMULATED® 52 


WEEK 43 


MO. OF ARRIVALS DURING PERIOD ® 9193 

MO. OF MSGS COMPLETED THIS PD.® 9192 
AVERAGE TIME IN SYSTEM ■ 1.6 

AVERAGE NO. IN SYSTEM ® 2.4 

AVERAGE FRACTIONAL CHAN. UTIL.® ... . 0.14 
AVERAGE TIME IN QUEUE ® 0.0 

AVERAGE PROCESSING TIME , ® 1.6 


THE MAXIMUM DELAY OF 0 MINUTES OCCURRED 9193 TIMES 




48 


REFERENCES 

1. Hillier, Fredericks.; and Lieberman, Gerald J. : Introduction to 

Operations Research. Holden-Day, Inc., 1967. 

2. Miller, Irwin; and Freund, JohnE.: Probability and Statistics for 

Engineers. Prentice -Hall, Inc. , 1965. 

3. Carnahan, Brice; Luther, H. A.; and Wilkes, James O. : Applied 

Numerical Methods. John Wiley & Sons, Inc., 1969. 



TABLE I. - TABLE OF SEASONAL VARIATIONS OF DISASTER WARNING MESSAGES FROM 1966 - 1971 











































FOR CHANNELS NUMBERING 1 TO 20 


Number 

channels 

Average 
time in 
system 

Average 
number 
in system 

Average 

fraction 

channel 

utilization 

Average 
time in 
queue 

Average 

processing 

time 

Maximum 

delay 

1 

2.9 

3.7 

0. 77 

1.4 

1. 6 

13 

2 

2.9 

3.7 

0.74 

1.4 

1. 6 

13 

3 

1.8 

2. 6 

0.49 

0.2 

1.6 

5 

4 

1.6 

2.4 

.37 

.0 

1. 6 

4 

5 

1.6 

2.4 

.30 

.0 

1. 6 

2 

6 

1.6 

2.4 

.25 

.0 

1.6 

1 

7 

1.6 

2.4 

.21 

.0 

1.6 

1 

8 

1. 6 

2.4 

. 19 

.0 

1. 6 

1 

9 

1.6 

2.4 

. 16 

.0 

1. 6 

0 

10 

1.6 

2.4 

.15 

.0 

1.6 

0 

11 

1.6 

2.4 

. 13 

.0 

1.6 

0 

12 

1.6 

2.4 

.12 

.0 

1.6 

0 

13 

1.6 

2.4 

.11 

.0 

1.6 

0 

14 

1.6 

2.4 

.11 

.0 

1. 6 

0 

15 

1.6 

2.4 

.10 

.0 

1. 6 

0 

16 

1.6 

2.4 

.09 

.0 


0 

17 

1.6 

2.4 

.09 

.0 


0 

18 

1.6 

2.4 

.08 

.0 

1. 6 

0 

19 

1.6 

2.4 

.08 

.0 

1. 6 

0 

20 

X.6 

2.4 

.07 

.0 

1. 6 

0 


Random number seed = 8 

Number of arrivals = 9357 

Number of messages completed = 9355 












TABLE IV. - PROBABILITIES OF BEING 


IN STATE ZERO AND STATE (S+l) 
FOR X = 0.9717 AND ju = 1.008 


s 

P 

0 

P S + 1 

3 

0.377555668 

0.026689434 

4 

.380904197 

.004351790 

5 

.381315291 

.000632013 

6 

.381363153 

.000081526 


.381368398 

.000009239 


.381368994 

.000000623 


.381369114 

<.0000001 

10 

.381369114 

<. 0000001 

11 

.381369114 

<. 0000001 

12 

.381369114 

<. 0000001 

13 

.381369114 

<. 0000001 

14 

.381369114 

<.0000001 

15 

.381369114 

<. 0000001 

16 

.381369114 

<.0000001 

17 

.381369114 

<.0000001 

18 

.381369114 

<* 0000001 

19 

.381369114 

<.0000001 

20 

.381369114 

<.0000001 










TABLE V. - PROBABILITIES OF BEING 
IN STATE ZERO AND STATE (S+l) 


FOR A = 0.9717 AND n = 0.625 


s 

P 

0 

p s*l 

3 

0. 197496175 

0.106376171 

4 

.208861827 

.027976811 

5 

.210840284 

.006570279 

6 

.211182415 

.001367569 

7 

.211238622 

.000253737 

8 

.211247265 

.000042796 

H 

.211248517 

.000006914 

m 

.211248695 

.000001430 

M 

.211248755 

<. 0000001 

12 

.211248755 

<.0000001 

13 

.211248755 

< 4 0000001 

14 

.211248755 

<; 0000001 

15 

.211248755 

<.0000001 

16 

.211248755 

<. 0000001 

17 

.211248755 

<. 0000001 

18 

.211248755 

<.0000001 

19 

.211248755 

<. 0000001 

20 

.211248755 

<. 0000001 









TABLE VI. - PROBABILITIES P g+1 FOR VARIOUS 
PERCENTAGES OF DEGRADATION FOR 


S = 10, X = 0.9717, (i = 1. 008 ; 
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